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ABSTRACT 

Future observations with ALMA will be able to detect star-forming primordial galax- 
ies, and perhaps even their dwarf spheroidal satellites, in CO emission lines at redshifts 
approaching 10. These observations will compliment other tools designed to study the 
dawn of galaxy formation, such as NGST and FIRST. 
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1 INTRODUCTION 

Modern cosmology continues to push the boundaries of the 
known universe to higher and higher redshifts. We are de- 
signing telescopes capable of detecting the very first objects 
that formed in the universe. The real breakthrough will oc- 
cur in the next decade: the New Generation Space Telescope 
(NGST), the The Far IR and Submillimiter Space Telescope 
(FIRST), and the Atacama Large Millimeter Array (ALMA) 
will allow us to reach cosmological redshifts in excess of 10. 
ALMA is especially important, because it can observe CO 
emission from primordial galaxies redshifted into the mil- 
limeter band. Because the Cosmic Microwave Background 
(CMB) temperature increases toward the past, higher ro- 
tational l evels of the CO mo l ecule are popul a ted at high 



redsh ifts (Silk fc Spaans 1997 
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§00$, 

resulting in a large negative K-correction. The effect 
is so large that a star-forming galaxy will appear equally 
bright at z — 5 and at z = 10. 
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Frayer et al. 1998; Downes et al. 1999|; |Frayer et al. 


1999 


Andrcani et al. 200C; Papadopoulos et al. 2000|; |Pa- 


padopoulos et al. 2001). These results provide information 



on the kinematics and energetics of the star-forming (molec- 
ular) ISM and thereby complement observations of the stel- 
lar light component in primordial systems. 

In this paper we show that ALMA, with its significantly 
higher sensitivity, will be able to image primordial galaxies 
to redshift 10 and even beyond. It is conceivable that it 
will even be able to detect dwarf spheroidal satellites of pri- 
mordial galaxies and their tidal tails (which should be in 
abundance at high redshift, where the merger rate is much 
higher than in the local universe). 



Table 1. Simulation Parameters 

Mass res. Spatial res. 
10 1.0/i _1 kpc 



Run N 
A 128 3 
B 64 3 



Box size 
4h _1 Mpc 
2/1" 1 Mpc 



10 6 - 6 M r . 
10 6 - 6 M, 







1.5/i _1 kpc 



2 METHOD 
2.1 Simulations 

We use the cosm ologic al simulations of reionization re- 
ported in Gnedin ( 2000| ) . The simulations include 3D radia- 
tive transfer (in an approximate implementation) and other 
physical ingredients required for modeling the process of cos- 
mological reionization. 

Two simulations of a representative CDM+A cosmolog- 
ical model^jwere performed with the parameters specified in 
Table |l| Both simulations were stopped at z = 4 because at 
this time the rms density fluctuation in the computational 
box is about 0.25, and at later times the box ceases to be a 
representative region of the universe. 

The two simulations from Table ^ allow us to investigate 
the sensitivity of our results to the inevitably missing small- 
scale and large-scale power. The difference between the two 
runs can be interpreted as the theoretical uncertainty in our 
calculations, given a cosmological model. Clearly, our results 
will be different for different assumptions about cosmological 
parameters. 



2.2 Calculating CO emission from early galaxies 



The code d escrib ed in Spaans (1996) and applied as in Silk 



Spaans (1997) has been used to rerun the models pre- 



With the following cosmological parameters: f2o = 0.3, fl\ = 
0.7, h = 0.7, fl b = 0.04, n = 1, a* = 0.91, where the amplitude 
and the slope of the primordial spectrum are fixed by the COBE 
and cluster-scale normalizations. 
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sented in ( |Silk fc Spaans 1997| ) with the latest atomic and 
molecular collision and chemistry data. These models use 
the Orion molecular cloud and its so-called bar region as be- 
ing representative of a region of active star formation. The 
star formation rate in Mq yr" 1 of a fiducial model galaxy is 
then related to the total number of Orion-like star formation 
sites through division by the average star formation rate of 



the Orion region, ~3x 10~ 4 M© yr" 1 (Hillenbrand 1997) 



In Silk & Spaans (1997) it has been shown that the CMB 



becomes an important source of excitation at high redshift 
because of the fortuitous coincidence between the CO level 
spacing and the 1 + z increase in the CMB temperature. 
This causes galaxies at z = 5 and z = 10 to be observable at 
similar flux density levels, provided they in fact are present. 
It has been assumed that the Orion-like regions responsible 
for the star formation activity occur throughout the model 
galaxy, and are not all conf ined t o the inner few 100 pc as in 
Combes, Maoli, & Omont (1999) This assumption decreases 
the mean optical depth of the CO lines and is most likely to 
hold at high [z > 3) redshifts, when galaxies are still being 
built up through mergers and accretion. 

In order to compute the spectrum of CO emission as a 
function of wavelength, for a given bandwidth Ai < A < A2 
and a given transition J — » J — 1, we identify a range of cos- 
mological scale factors ai < a < ai that correspond to our 
bandwidth. This range of scale factors in turn corresponds to 
the range of comoving distances x\ < x < X2- However, due 
to periodic boundary conditions adopted in the simulations, 
we cannot always model this range of comoving distances 
directly - if it is large enough, it will correspond to more 
than one box size. If we simply stack a sufficient number 
of simulation boxes together, we will get an unphysical re- 
sult due to periodicity. In order to break this perio dicity , 
we use the approach described in Gnedin & Jaffe (2001): 



we randomize the neighboring boxes by randomly flipping, 
transposing, and shifting each of the periodic images of the 
computational box. 



3 RESULTS 

Figure |l| serves to illustrate the uncertainty of our calcu- 
lations due to the finite size of the computational box and 
finite numerical resolution, as measured by the difference be- 
tween the two simulations A and B. In addition, since both 
simulations were stopped at z = 4, a contribution from later 
redshifts cannot be included. In order to estimate the effect 
of this contribution, we calculated the CO emission for two 
cases: no star formation after z = 4, and constant star for- 
mation after z = 4. The difference between those two cases 
quantifies the uncertainty due to the finite value for the final 
redshift of our simulation. 

For A < 0.3 cm our calculation is not reliable even 
in a qualitative sense (to within a factor of 2). At higher 
wavelengths finite numerical resolution still prevents us from 
achieving better than about 50% accuracy. More than that, 
since the star formation rate in our simulations is normalized 
to the observed value at z = 4, which is in turn uncertain to 



at least a factor o f two ( Nagamine, Cen, fc Ostriker 2000 
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Figure 1. The mean CO flux density as a function of wavelength 
from two simulations: the small simulation box (dashed line) and 
the large box (solid line) , as described in Table 1 . Thin lines show 
the calculation with the star formation rate set to zero after the 
final redshift of simulation (2 = 4); bold lines are for the case 
when the star formation rate is assumed constant for z < 4. 
The difference between the two sets of lines is our theoretical 
uncertainty due to the final size of the simulation box. 



than to make some definite predictions. Because ALMA will 
not be commissioned until about 2010, theorists have plenty 
of time to improve upon our calculations and come up with 
more definitive predictions. 

Figure ^ shows the evolution of the mean mass- weighted 
metallicity of the gas and stars in our large simulation (run 
A). One can see that stars quickly gain a metallicity of 3% 
solar by z ~ 15, and then gain another order of magnitude 
on average by the end of the simulation at z = 4. The metal- 
licity of the gas is always lower than stellar, but increases 
more rapidly. The decrease in stellar metallicities at higher 
redshifts slightly mitigates the increase in the CO emission 
due to higher CMB temperature, however it is not sufficient 
to completely remove the negative K-correction, and thus 
high redshift star forming galaxies should be considerably 
brighter than their low redshift counterparts. 

In Figure ^ we show our main result: the 4 square de- 
grees of the sky (an image of our computational box) at 
four different wavelengths. The peak of CO emission corre- 
sponds to a broad r ange of transitions w ith values in the 
range of J = 5 to 12 (Silk & Spaans 1997). This means that 



Steidel et al. 1999), our results in general are uncertain to 
a factor of two to three. However this is quite sufficient for 
our purpose, which is to emphasize the possibilities rather 



several epochs correspond to a single observed wavelength: 
A = 0.3 cm maps redshifts from about 5 to 13, A = 0.6 cm 
corresponds to redshifts from about 10 to 30, and A = 1.0 cm 
includes everything from z — 20 to z = 50. This is some- 
what unfortunate, because it means that a simple image at 
a given wavelength can only provide information about star 
formation over a range of redshifts. On the other hand, pri- 
mordial star forming galaxies have relatively narrow velocity 
ranges, which means that by observing the CO emission in 
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Figure 2. Evolution of the mean mass-weighted metallicity of 
the gas (dashed line) and stars (solid line). 
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Figure 3. A 2 X 2 degree image of the sky at four different wave- 
length (A = 0.3, 0.4, 0.6, and 1 cm clockwise from the top left 
panel, as labeled in the figure) taken with a frequency resolution 
R = lg(A/AA) = 3. Side bars give the correspondence between 
the color and decimal log of flux density in mjy per square arc- 
second. 



a narrow wavelength band, it is still possible to disentangle 
different epochs. 

This is illustrated in Figure ^, where we show the same 
wavelength but at four different spectral resolutions. With 
a resolution A/ A A = 10 4 we can eliminate most of the emis- 
sion from different redshifts and obtain an image of a single 
primordial galaxy (lower left panel in Fig. ^|). Notice that 
not only the galaxy itself is measurable by ALMA, but also 
its satellites - young dwarf spheroidals - are clearly visible 
in the image. A fortunate circumstance - an enhancement in 
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Figure 4. A 2 X 2 degree image of the sky at A = 0.3 cm taken 
with four different wavelength resolutions R = lg(A/AA) = 1, 
2, 3, and 4, clockwise from the top-left panel (as labeled in the 
figure). Side bars give the correspondence between the color and 
decimal log of flux density in mjy per square arcsecond. 
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Figure 5. Fraction of the sky above a given flux density for 
three different spectral resolutions at A = 0.3 cm. The black ar- 
row shows ALMA's sensitivity per 1 arcsec 2 pixel for a 100 hour 
integration (or, equivalently, 1 hour integration on a 10 arcsec 2 
source) . 



population of high J levels due to higher CMB temperature - 
may allow ALMA even to see dwarf spheroidals being tidally 
disrupted in the vicinity of larger galaxies and to map their 
tidal tails (although this will be a very difficult observation). 
And, of course, the high redshift picture is much more spec- 
tacular than its local equivalents because the merging rate 
at z ~ 10 is some 500 times higher than in the local universe. 
The distribution of flux density on the sky is shown in 
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Figure 6. Fraction of the sky above the ALMA sensitivity as a 
function of the wavelength. We only show wavelengths A > 0.2 cm 
since at shorter wavelengths our calculations become unreliable, 
as shown by Fig. [l]. 

Figure ^. We also show a sensitivity limit for ALMA for a 
typical observation. While ALMA will be just short of mea- 
suring typical objects with a spectral resolution of 10 3 , at 
a resolution of 10 4 it will be able not only to detect pri- 
mordial galaxies but also to take their images. High spectral 
resolution is also required to separate emission from differ- 
ent redshifts, so this mode appears to be the most promising 
for high redshift observations with ALMA. 

Figure ^summarizes our results by giving the fraction of 
the sky above the ALMA sensitivity as a function of wave- 
length. ALMA will have an easy time observing primeval 
star-forming galaxies at z ~ 5 — 10, but will not be able 
to see much beyond that because of the strong decrease in 
stellar metallicities and star formation rates at higher red- 
shifts. The latter statement is, of course, strongly dependent 
on the cosmological model - models with large amounts of 
small scale power will have star formation commencing ear- 
lier. 



4 CONCLUSIONS 

We have demonstrated that future observations with ALMA 
will be able to detect star-forming primordial galaxies, and 
even their dwarf spheroidal satellites, in CO emission lines, 
mostly due to a large negative K-correction. High spectral 
resolution observations are required to both separate contri- 
butions from several objects at different redshifts all emit- 
ting at the same wavelength from different rotational levels, 
and to increase the signal-to-noise. Direct imaging will be 
possible for galaxies up to redshifts approaching 10. Unfor- 
tunately, the expected decrease in the star formation rates 
and stellar metallicities at higher redshifts will make obser- 
vations of the era before z = 10 extremely difficult. 



A simple possible strategy for these observations might 
include a shallow large area survey to identify early forming 
massive galaxies (rare peaks), with subsequent targeted deep 
observations of their environments. In this case it is possible 
to observe the same region on the sky at several wavelength 
bands spaced around several subsequent CO emission lines, 
in order to get a linear increase in the signal-to-noise. 

Broadly speaking, CO emission lines, boosted by higher 
CMB temperature in the dawn of galaxy formation, will al- 
low millimeter observations with ALMA and other instru- 
ments not only to complement observations with the Next 
Generation Space Telescope (NGST), but, perhaps, to di- 
rectly compete with the NGST in the race for the highest 
observable redshifts. 

This work was partially supported by National Com- 
putational Science Alliance under grant AST-960015N and 
utilized the SGI/CRAY Origin 2000 array at the National 
Center for Supercomputing Applications (NCSA). 
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